Quantum Dot in 2D Topological Insulator: The Two-channel Kondo Fixed Point 
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In this work, a quantum dot couples to two helical edge states of a 2D topological insulator 
through weak tunnelings is studied. We show that if the electron interactions on the edge states are 
repulsive, with Luttinger liquid parameter K < 1, the system flows to a stable two-channel fixed 
point at low temperatures. This is in contrast to the case of a quantum dot couples to two Luttinger 
liquid leads. In the latter case, a strong electron-electron repulsion is needed, with K < 1/2, to reach 
the two-channel fixed point. This two-channel fixed point is described by a boundary Sine-Gordon 
Hamiltonian with a K dependent boundary term. The impurity entropy at zero temperature is 
shown to be In v2_R'. The impurity specific heat is C ccT i< when 2/3 < A' < 1, and C ccT when 
K < 2/3. We also show that the linear conductance across the two helical edges has non-trivial 
temperature dependence as a result of the renormalization group flow. 



Introduction — The two-channel Kondo model has been 
under intense theoretical studj- ^ after the seminal work 
of Noziere and Blandin'^ , in which they pointed out that 
a non-Fermi liquid fixed-point exists for the two-channel 
Kondo model with non-interacting electrons. Recently, 
the authors of Ref 3 demonstrated experimentally that a 
quantum dot couples to an infinite reservoir with non- 
interacting electrons and a finite reservoir with interact- 
ing electrons may display two-channel Kondo effect. The 
key idea is that Coulomb blockade suppresses the ex- 
change of electrons between the finite reservoir and the 
infinite reservoir.!^ Consequently, the two reservoirs cou- 
ple to the quantum dot independently. The two-channel 
Kondo quantum critical point can be reached by tuning 
the couplings of the quantum dot with the reservoirs to 
equal strength. 

Indeed, using electron-electron repulsion to suppress 
interreservoir tunnelings was suggested earlier by Fab- 
rizio and Gogolin, in the context of a quantum dot cou- 
ples to two Luttinger liquid leads.l^ They show that 
the two-channel Kondo fixed point can be reached when 
the repulsive interactions of the electrons in the leads 
are strong enough, i.e. with Luttinger liquid parame- 
ter K < 1/2. Similar conclusions were also made in later 
works z^ *" However, two-channel Kondo effect in Luttinger 
liquids have never been observed experimentally. 

Recently, a new class of materials called topological in- 
sulators were first th eoretically proposed and then exper- 
imentally fabricated.'^'^^^^^^^^Two dimensional topolog- 
ical insulators have gapless helical edge states despite the 
presence of a bulk gap .1^ The edge states are called hehcal 
because the directions of the spins and the momentum of 
the electrons are tied together.— The one-channel Kondo 
effect of helical edge states was first studied in Ref ITU as 
a truncated model of Kondo effect in Luttinger liquids. 
A more detailed recent study is carried out in Ref 15. 

In this work, we study a quantum dot couples to two 
helical edge states as depicted in Fig. la. If the electrons 
on the edge states are non-interacting, it is well known 
that the syst em ca n be described by a one-channel Kondo 
Hamiltonian.'^^'^^' However, we show that a weak repul- 
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FIG. 1: a) A quantum dot couples to two helical edges of 
a 2D topological insulator through electron tunneling. S and 
D denote the source and drain respectively. b)The schematic 
picture of the linear conductance G from SI to D2 versus 
temperature T. The solid line and the dashed line depict the 
1/2 < A < 1 and K < 1/2 cases respectively. 



sive interaction, with if < 1 is enough to drive the system 
to the two-channel Kondo fixed point. This is in sharp 
contrast to the case of Luttinger liquid leads in which 
K < 1/2 is needed.!^ This two-channel fixed point is de- 
scribed by a boundary Sine-Gordon Hamiltonian with a 
K dependent boundary term. We study the ground state, 
thermaldynamic and transport properties near this two- 
channel fixed point. More specifically, we s how that the 
impurity entropy at zero temperature is In \J1K, the im- 
purity specific heat Cimp oc T#~^ when 2/3 < K <\, 
and Cirap oc T when K < 2/3, and the linear conduc- 
tance across the two helical edges through the quantum 
dot has non-trivial temperature dependence as a result 
of the renormalization group fiow. 

The Model — The schematic diagram of our model is 
depicted in Fig. la. Two helical edge states of an 2D 
topological insulator are brought close to each other at 
a tunneling junction. A quantum dot is placed at the 
middle of the junction and couples to the two edge states 
through electron tunnelings. The system is described by 
an Anderson Hamiltonian Ha = Hq + Hi + Hd + Ht where 
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and Hi denote the kinetic energy and the interac- 
tion of the edge states respectively. Hd is the Hamilto- 
nian of the quantum dot where e is the single particle 
energy level at the dot and U is the on-site interaction 
energy. Ht describes the couplings between the dot and 
the edges and we have tuned to the channel symmetric 
point such that t is independent of the edge index. In the 
presence of spin orbit coupling, spin ("f and J,) denotes 
pseudo-spin of the edge modes. In the large U limit, one 
obtains the Kondo Hamiltonian H — Hq + Hk through 
perturbation theory,^® where 



4=1,2 i^ftj 

. . . (1) 

Initially, Ji = J2 > and without intera ction , Eqjl] re- 
duces to a single channel Kondo probleml ^^ l ^^ l Our goal 
is to show that J2 normalizes to zero at low temperatures 
when 52 > 0. Because of helicity, one of the R/L and the 
spin "f I indexes are redundant, we keep only one of them 
starting from Eq[T| 

It is important to note that, in general, there are ex- 
tra backscattering terms of the form ?/'|fl(O)'0jL(O) -I- h.c. 
caused by the quantum dot in the Hamiltonian H. The 
backscattering terms cut each edge into two separate 
parts at low temperatures if if < 1 and render the ge- 
ometry in Fig. la unstable. However, in our case, sin- 
gle particle back-scattering terms are forbidden because 
of time-reversal symmetry. Two-particle backscatterings 
preserve time-reversal symmetry but they are irrelevant 
so long as if < 1/4.^"^ Thus, our discussion below is valid 
for if > 1/4. 

By Abelian Bosonization ,112120] the electron operators 

1 p±i{\/4TTcf>ip,j[^(x)+kFx) 



can be written as: ipiR/L = - 
We define Bosonic fields (j>i{x) = 4>iL{x) -\- 4>iR{x) and 
6*^(2;) = (piLix) — (j)iu{x). The dual fields satisfy the 
commutation relations [(j)i{x),9j(x')] = ^SijSgii(x — 
x'), where sgn(a;) = when a; = 0. The sym- 
metric and antisymmetric combinations of (j)^ and 9i 
are denoted as: (pg/^ = ~^{4'i i 4'2) and Og/^^ = 
± 62). In terms of (ps/a and Og/a, the 
above Hamiltonians can be written as: Hq + Hi — 
^ / j;{d,cPg)^+K{d:,9g)^ + j^{d^cba)^+K{d,0a)^dx and 
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Here, 



HK = -^^iJ!S,d^0g{O) 

■ '^'""'-(S-e-'^"*=(o) -KS'+e'^*=(o))cos(\/2^0a(O)) 
-5*^ sin(V2^6'a(0)) sin(^/2^0a(O)) 

-g-zV2^0,(O) ^ s+^^V2^<^A0)) cOs(\/2^0a(O)). 

the Luttinger liquid parameter is de- 
K = 
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as 
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VF^/i^ + g4/2TrvF)^ - (32 /27rt)p-) 2.1221 The Bosonic 
form of Hii is spin anisotropic, we denote the two com- 
ponents of Ji as J^^^^ respectively. This is reasonable 
for our problem, because spin-orbit coupling breaks 
SU{2) symmetry of the electrons on the edges. The 
spin-coupling terms acquire different scaling dimensions 
when if ^ 1 as we see below. 

Renormalization Group Analysis — With the Bosonic 
Hamiltonian H = Hq + iij + H^ , we may calculate the 
scaling dimensions of the operators in Hk- At the vicin- 
ity of the fixed point where Jf = ^ = 0, the scaling 
dimensions of the Jf , Jf^ terms are 1 and if respectively. 
On the other hand, the scaling dimensions of the J| and 
terms are 1{K + 1/if) > 1. Thus, when if < 1, the 
J2 terms decrease but the Ji terms grow when tempera- 
ture is lowered. 

In order to study the physics at the strong cou- 
pling regime when Ji is of order 1, we use the Emery- 
Kivelson method by applying a unitary transformation 
U = e'v^2^*'(o)'S" to iJ.ISil After the unitary transforma- 
tion and the rescaling of the bosonic fields, we have: 

ii" = Whu 

= Ho + XS^d^dgiO) + ^S^sm{^0amMV2^MO)) 
+ S+)l^ cos(y2^0,(O)) + ^ cos(y^0,(O))], 

(2) 

where A = \/2iiKv'p — \J^^Ji- At the vicinity of the 

point A = J^^ = J| = J^^ = 0, we may calculate the 
scaling dimensions of the Ji terms as before. The J^^ and 



tion of the e**'^'''-'^' factors by the transformation. The 
J^^ term has scaling dimension ^. It is relevant when 
if < 1. The J 2^ term has scaling dimension and it 



is relevant when 1/2 < if < 1. 

If (pa and 9a are independent fields, both J^^ and J, 
flow to infinite. As a result, the two-channel fixed point 
cannot be reached when 1/2 < if < 1. That is the 
case for a quantum dot couples to two Luttinger liquid 
leads.^ However, in the present problem, (pa and 9a are 
dual fields of each other. This fact changes the physics 
quite dramatically. 

When if < 1, the Jf^ term is the most relevant term. 
At low temperatures, one expect the effective Hamilto- 
nian of H is described by the fixed point Hamiltonian: 

2 /^^ 

Heff = H„ + ^S'=cos{V2T^M0))- (3) 
ira 

To show that this fixed point is indeed stable, 
we need to calculate the scaling dimensions of 
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the operators J^^ S"" cos{^J^ea{Q)), 5XSA0s{Q) and 

J|S';,sin(y^6'a(O))sin(V2^0a(O)) at the vicinity of 
the fixed point. 5\ denotes the deviation from the 
Emery-Kivelson hne. The results are given in the third 
column in Table I. We see that all these operators are 
irrelevant if < 1. Which confirms that Eq|3]is the low 
temperature effective Hamiltonian. Below we give some 
details of the calculation. 

To calculate the scaling dimension of 
Jj^S^ cos (^6*0(0)), we observe that commutes 
with Heff, so that we may set to ±1/2. As a 
result, one needs to calculate the scaling dimension of 

cos(a /^6'a(0)) with the Hamiltonian 



H±=Ho± 



J 



xy 



COS 



(72^0,(0)). 



(4) 



Eq|4]has the form of the Hamiltonian of a spinless Lut- 
tinger liquid wire with an impurity backscattering term 
at point X — 0. One may interpret the cos(\/27r^0a(O)) 
as an impurity backscattering term with Luttinger liq- 
uid parameter K' = K/2. In a spinless Luttinger liq- 
uid, the backscattering term is relevant when K' < 
1. At low temperatures, the "backscattering term" 
cos{\/27rK(j)a{0)) cuts the Luttinger liquid wire into two 
separate pieces at a; = 0.1^ As a result, the point a; = is 
not a point at the bulk anymore, it can be regarded as a 
point at the boundary when J^^ ^ 00. In the bulk, the 

scaling dimension of cos \J^&a{0) is 1/2K. However, at 

the boundary, 0q(O) is pinned to a constant value when 
J^^ — !■ 00. Thus, the value of (t>a{x) for a; — > is con- 
strained to be near the value of 0a(O), by the {dx<pa{x))'^ 
term in Hq. Hence, the fluctuation of 0^(0) is enhanced. 
Detail calculations show that the scaling dimension of 

cos ^J^6a{0) should be 1/K instead.'^S! Consequently, 

the cos ^ ^Oa{0) term is irrelevant when K < 1. 

The calculations of the scaling dimension of the 
S^dxOs(fi) term is similar to what have been done before 
in Ref|23] Since 9s is governed by a free Hamiltonian, we 
have < S''{T)dx9s{T)S^{0)dxes{0) >=< S'^(r)S'^(0) >< 
dx0s{T)dx9M > 

In order to calculate < S'^(r) 5*^(0) >^ , we note that 
is the sum of the raising and lowering operators of the 
eigenstates of S"^, |± >. Thus, we have 



< -|e^'=//^S'^(0)e--^''//^S'^(0)|- > 
=< -|e^-^M/^t(o)e-^-^T4^(0)|- > 
^<-\WHT)Wm->f,_ 



(5) 



where W is a unitary transformation which satisfies 



w 



-H+T 



Actually, we need to find a unitary operator W 
which transforms 0^ to 0^ + ^ Indeed, W — 
55f (ea(L)-e„(-L))^ where 2L is the length of the sys- 
In order to calculate < — |W^'f (t)VF(0)|— >^_, we 
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TABLE I: Scaling dimensions of the operators (Op.) at dif- 
ferent fixed points (FPs). (J2 = for all fixed points, only 
coefficients of corresponding operators sliown.) 



Op. \ FPs 


Ji = 


A = = 


A = 0, 00 


Jf(5A) 


1 


1 


1 + 1/2K 




K 


K/2 


N/A 


J2 


UK + 1/K) 


UK + 1/K) 


1/K+1/2K 




UK + 1/K) 


1/2K 


1/K 


V2 


UK + 1/K) 


UK + 1/K) 


1/K 



come back to the spinless Luttinger liquid Hamiltonian 
interpretation of 

Suppose we have a spinless Luttinger liquid wire with 
an impurity at a; = and the left and right ends lo- 
cate at — L and L respectively. At the strong coupling 
fixed point, the wire is cut into two separate parts by the 
impurity. However, the left half and the right half are de- 
scribed by free Hamiltonians Hql and Hqh respectively. 
As a result, we have 



< W^t)W{0) >h 

(0„(-L,r)-e(-L,O)) 



cx 



< e' 
1 



.(e„(L,r)-e(L,0)) 



(6) 

Since dx9s{Q) has a scaling dimension of 1, we conclude 
that the operator S^dx9s{0) is irrelevant with scaling di- 
mension 1 / 2K + 1 . 

One may calculate the scaling dimension of the 
JfS'z sm{^J^9a{0)) sin(-\/27rif</)o(0)) term in a similar 

way. This term has scaling dimension j]^ + It is 
irrelevant when K < I. 

Up to now, we have neglected the J2 V^t/'l^ipja term 

in Hk- However, this term is always irrelevant and we 
will drop this term for the rest of the paper. 

Impurity Entropy — We have shown above that the 
two-channel fixed point is stable. At the vicinity of this 
fixed point, the system is described by Heff in Eqjsj It 
is interesting to note that H^ff can be regarded as a 
generalization of the fixed point Hamiltonian of the two- 
channel and four-channel Kondo non-Fermi liquid fixed 
point Hamiltonians. The two-channel and four-channel 
cases are described by H^ff with K — 1 and if = 3/2 
respectively.'^^ One of the most remarkable features of the 
multichannel Kondo effect is the existence of fractionally 
degenerate ground state. For example, the two-channel 
and four-channel Kondo mod els ha ve residual entropy of 
I In 2 and ^ In 3 respectively.'^^'^ This motivates us to 
study the residual entropy of Hf.ff for if < 1. 

In order to calculate the residual entropy, we first 
calculate the partition function Z = Trj exp^'^^^-f^ }. 

Since commutes with the H^ff, the states can be la- 
beled as I ± 



>HaE 



5'" 



>. We may write Z = Z+ + Z. 



4 



Tr{exp-'^^+ } + Trjexp"^-'^- }. Remember that the 
integral of Hq run from — oo to +00 if L ^ 00 is taken. 
However, we may define new bosonic fields such that 
the positive axis is folded to the negative axis!^ After 
folding, the effect of the cosine term in H± is to intro- 
duce a Dirichlet boundary condition. The actual value 
of (f>a{0,T), whether it is pinned at {2n + l)n/\/2TTK 
or 2mr/y/2TTK, does not affect the partition function.'^® 
Thus, we have = Z_ and Z = 2Z^. The impu- 
rity entropy of has been calculated in Refl27l It is 
In g = In y^K/2, where g is the "ground state degener- 
acy" discussed in RefUS) Together with the contribution 
from , the total impurity entropy is 

= In V2K. (7) 

Substituting K = 1 and K = 3/2 into Eqj?] repro- 
duces the impurity entropy results, calculated by the 
Bethe ansatz and boundary conformal field theory, for 
the two-channel and four-channel Kondo effects with 
non-interacting electrons respectively .'^*^ 1^ 

Specific Heat — We show above that the residual impu- 
rity entropy is determined by the boundary term. On 
the other hand, the thermodynamic properties are de- 
termined by the leading irrelevant operators (LIO) near 
the fixed point. Since the irrelevant operators have 
correlation functions of the form << O(t)O(0) >>= 
l^^^^i^]^^, the second order correction of the free en- 

Lsin(7rTT)-l ' 

ergy SPi^^T) = -A^ J^J^ dri^^^^]^^ where tq is the 

cut-off time.f^ 

Near the strong coupling fixed point, the 

cos{\/ ^6a{0)) term with scaling dimension 1/K is 

the LIO if 1/2<K < 1. We have 6F cx T^/^^-^ + 0{T^), 
where the T^/if-i ^gj.jjj jg universal and independent 
of the cut-off time tq and the O(T^) term has a tq 
dependent coefficient. As a result, the impurity specific 
heat C„ip cx T^/^-^ if 2/3 < K < 1. However, when 
K < 2/3, the free energy is dominated by the 0(T2) 
term and Cj^p oc T. 

Linear Conductance — In the above calculations, we 



have assumed that the two helical edges have equal chem- 
ical potential. If a small bias is applied across the two 
edges, there should be a current fiowing from one edge 
to the other. At temperatures much higher than the 
Kondo temperature, the linear conductance across the 
edges can be written as G cx ^J^^)^ + ( J|)2]r^+i/^f-2 
by pertubation to second order. At low temperatures, 
the conductance can be obtained by substituting the nor- 
malized value of J2 into the high temperature expres- 
sion of G. As we show above, J2 decreases near the 
fixed point with J^^ — ~ J2 ~ and Ji grows. If 
1/2 < < 1, however, when the Ji terms are of order 
1, the ^ term becomes relevant and grows. When the 
temperature decreases further, the strong coupling fixed 
point with J^^ — s- cx3 is reached. Near this fixed point, 
J2 ^ and J| are irrelevant and fiow to zero. As a result, 
the conductance acquires a non-trivial temperature de- 
pendence. On the other hand, if if < 1/2, the J2 terms 
are always irrelevant and decrease to zero. Consequently, 
the conductance decreases monotonically as the temper- 
ature lowers. The temperature dependence of the linear 
conductance at different fixed points are shown in the 
Fig.lb. 

Conclusion — We show above that a quantum dot cou- 
ples to two helical edge states of 2D topological insulators 
is described by a Kondo Hamiltonian. Weak repulsive in- 
teraction with if < 1 on the edges, drives the system to a 
two-channel fixed point at low temperatures. At the fixed 
point, the residual impurity entropy is In \J2K. Near the 
fixed point, the impurity specific heat is Cimp cx r2/if-2 
when 2/3 < if < 1 and C,„p cx T when K < 2/3. 
Moreover, the linear conductance across the two edges 
through the quantum dot has non-trivial temperature de- 
pendence. 
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